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We study (2+l)-dimensional multicomponent spatial vector solitons with a nontrivial topological 
structure of their constituents, and demonstrate that these solitary waves exhibit a symmetry- 
breaking instability provided their total topological charge is nonzero. We describe a novel type of 
stable multicomponent dipole-mode solitons with intriguing swinging dynamics. 



Recent progress in the study of spatial optical soli- 
tons and their interaction, as well as the extensive exper- 
imental demonstrations of stable self-focussing of light in 
different types of nonlinear bulk media, open the road 
for new concepts to control the diffraction of optical 
beams and to design new devices for optical switching 
and storage Many novel fundamental concepts re- 
cently suggested in the physics of spatial optical solitons 
are associated with vectorial interaction and multicom- 
ponent soliton beams that mutually self-trap in a non- 
linear medium. Such composite multimode solitons can 
have complex structures and, in many cases, their total 
intensity profile exhibits multiple humps ||] . 

In a bulk medium, vector solitons exist in different 
forms and, as was recently shown for two-component self- 
trapped beams, many types of multipole vector solitons 
can be predicted and analyzed for an isotropic nonlin- 
ear bulk medium with saturable nonlinearity jsj. Re- 
cently, an important generalization of this concept to the 
case of A^— component two-dimensional vector solitons 
was suggested for an example of a thresholding nonlinear- 
ity Q. In particular, MussHmani et al. Q predicted the 
existence of multihump A'^— component composite spatial 
solitons that carry different topological charges ('spins') 
and, therefore, can provide exciting possibilities for 'spin- 
dependant' interactions of self-trapped optical beams [^. 

The purpose of this Letter is twofold. First, we study 
in more detail the dynamics of multicomponent spatial 
solitons carrying topological charges in different compo- 
nents and demonstrate that, in contrast to the conjecture 
of their stability made in Ref. [Q, these vector solitons 
demonstrate a symmetry-breaking instability in all the 
cases where their total angular momentum is nonzero. 
Second, based on earlier studies of two-component vec- 
tor solitons 1^ and the conceptual approach developed 
m Ref. |4|, we propose a novel type of stable multi- 
component vector solitons consisting of two perpendic- 
ular dipole components trapped by the soliton-induced 
waveguide. These vector solitons are studied here for 
the case of A^ = 3 components, which are shown to be 
the building blocks for the solitons composed of A^ in- 



coherently coupled dipole-mode beams ^ . Additionally, 
we demonstrate numerically that these novel vector soli- 
tons are very robust for a broad range of their parameter 
space, and they demonstrate intriguing swinging dynam- 
ics outside the stability domain, resembling long-lived ex- 
citations and vibrations of molecules. 

We consider the interaction of A^ mutually incoherent 
(2-t-l)-dimensional optical beams propagating in a bulk 
saturable medium, described by the normalized equa- 
tions [j = 1,2,..., AT), 
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where is the transverse Laplacian and z is the prop- 
agation coordinate. Equations (|^) describe, in a rather 
simplified isotropic approximation, screening spatial soli- 
tons in photorefractive materials 0| . 

To describe multicomponent vector solitons in the 
framework of the model first we look for stationary 



solutions in the form Ej(x,y,z) = Uj{x,y) exp{—iPjz), 
where Pj is the propagation constant and Uj{x,y) is the 
envelope of the j-th component. Then, introducing the 
dimensionless parameter Xj — (1 — /3j)/(l — and nor- 
malizing the field amplitudes, u 
coordinates, {x,y) — > {x,y)/^/Y 
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where / = P is the normalized total intensity, and 
F{I) = /(I -I- sl)~^ with the effective saturation param- 
eter s = 1 — /?!. 

First of all, following Musslimani et al. jj], we seek 
multicomponent radially symmetric solutions of Eqs. (j^) 
for which the main component ui{x,y) = Ui{r) has no 
nodes, but each of the components Uk {k > 1) carries a 
different topological charge, Uk{x,y) = Uk{r) eicp{imkO). 
We denote such states as (0, .., m^, ..), and an example 
for A^ = 3 is presented in Fig. 1(a), where the same in- 
tensity distribution corresponds to two different states, 
(0,+l,+l) and (0,-hl,-l). 

In order to study stability of these composite solitons, 
we propagate them numerically and find that, provided 
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the total angular momentum is nonzero, all these mul- 
ticomponent solitons undergo a symmetry-breaking in- 
stability and fragment into a number of the fundamental 
solitons, as shown in Fig. 1(b) for the case (0,-1-1,-1-1). 
This instability is similar to the instability of the vortex- 
mode solitons described earlier for the two-component 
model. The resulting incoherent superposition of two 
parallel dipole components U2 and U3 can be regarded as 
a generalization of a two-component dipole-mode soliton 
{ui, V} to a three-component solution {ui, 1*2,1*3} at 
-^2 ~ '^s with the help of a transformation of the dipole 
components, V {1*2,1*3}, where 1*2 — T^cosV' and 
M3 = V sinip (-0 is a transformation parameter). Such a 
straightforward generalization is indeed possible for the 
iV-component system (^). 
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FIG. 1. Evolution of the three-component soliton: (a) sta- 
tionary solution at 2 = 0, (b) the symmetry-breaking insta- 
bility of the (0, +1, +1) solution at z = 80 , (c) the long-lived 
quasi-stable propagation of the (0, -1-1, —1) state at z = 500. 
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FIG. 2. Families of the three-component dipole-mode soli- 
tons. (a) Soliton structure at A2 = 0.5 and A3 — 0.65, (b) the 
total and partial powers vs. A3 at fixed A2 = 0.5. 

The most important property of the (0, -f 1, — 1) solu- 
tion is that its total angular momentum is zero, and this 
makes it stable. In our calculations, this vector soliton 
was observed to be unchanged for the distances of the or- 
der of 10'^ diffraction lengths. However, being launched 
with additional noise, this soliton displays slowly growing 



modulations, as shown in Fig. 1(c). The total intensity 
of the modulated rings in Fig. 1(c) preserves the initial 
ring profile, resembling an incoherent superposition of 
two perpendicular dipole components |^ . While the vec- 
tor soliton, consisting of two crossed dipoles, has been 
shown to be unstable without the third main component 
, we found that the three-component dipole-mode soli- 
ton is stable in our numerical simulations. Stabilization 
of the vector ring in the presence of the third component 
can be explained by the physics of the soliton-induced 
waveguides. Indeed, two crossed dipoles, 1*2 and M3, rep- 
resent a vectorial guided mode of the induced waveguide. 
A nontrivial rotational transformation of such a solution 
(see Ref. |6) for details) allows to find a whole family 
of possible superpositions of these modes, including, as 
a particular case, the vortex components shown in Fig. 
1(a) and the iV-component dipole-mode soliton. 

In order to find the multicomponent solitary waves 
with a nontrivial geometry, we integrate the system (|^) 
numerically, by means of a relaxation technique, and find 
a novel class of the dipole-mode soliton that consists of 
perpendicularly oriented dipoles with different powers: 
the simplest possible solution of this type has = 3 
components, and it is described by two independent pa- 
rameters (A2,A3), as shown in Fig. 2(a). The family 
of these solitons ranges from solutions where the funda- 
mental mode dominates the entire structure to solutions 
where one of the dipoles dominates, as can be seen in 
Fig. 2(b), where for fixed A2 — 0.5 the power of the com- 
ponents Pj = / dv is shown as a function of A3. 

Numerical propagation of these solitons has shown that 
from the lower cutoff value for A3, where the intensity 
of the U3-component vanishes, up to a value of about 
A3 — 0.7 these vector solitons are stable, whereas for 
higher A3 they decompose to form different new struc- 
tures. As can be seen from Fig. 3, an unstable soliton 
breaks the symmetry along both symmetry axes of the 
initial distribution. The products of this instability (see 
the last row in Fig. 3) are a fundamental vector soliton 
and a rotating dipole-mode soliton, recently introduced 
in Ref. 1^ as 'a propeller soliton. Those two simpler 
solitons fiy away from each other after the break-up. 

Near the instability threshold, for 0.7 < A3 < 0.8, we 
observe very interesting and intriguing dynamics, associ- 
ated with weak oscillatory instabilities. Figure 4 shows 
a characteristic example of this dynamics, when the in- 
stability breaks the symmetry only along one of the sym- 
metry axes (parallel to the orientation of the stronger 
dipole). The product of this instability is a structure 
consisting of a tripole, a dipole and a nodeless beam. 
This structure is remarkably long-lived and it has, as the 
snapshots show, a swinging behavior resembling a swing- 
ing mode of a three-atom molecule. We could observe 
almost three periods of such oscillations, until a strong 
energy-exchange between the two dipole beamlets sets 
in and destroys this structure. We expect that the vi- 
brational degrees of freedom, which are likely associated 
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with long-lived soliton internal modes, should manifest 
themselves in the rich dynamics of soliton collisions, as 
it is known for the study of a two-component model H . 
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FIG. 3. Unstable propagation of a three-component 
dipole-mode soliton at A2 = 0.5 and A3 — 0.8, and its de- 
cay into a fundamental vector soliton and a vector propeller 
soliton (shown by dashed). 
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FIG. 4. Swinging dynamics of the vector dipole-mode soli- 
ton at A2 = 0.5 and A3 — 0.75. 

Having found these novel composite solitons for the 
isotropic model, we wonder if the similar multicompo- 
nent solitons can exist in an anisotropic nonlocal model 



which is more consistent with the experimentally stud- 
ied photorefractive nonlinearities ||lo[] . In order to ver- 
ify this, we have used the A^— component generalization 
of the Zozulya- Anderson model that takes into account 
the most important properties of photorefractive non- 
linearities and found similar classes of multicom- 
ponent localized solutions with perpendicularly oriented 
dipole components. Since the anisotropy allows stable 
stationary dipole modes oriented in two fixed directions 
only these solutions are found to be stable even in 
anisotropic media with nonlocal nonlinear response. This 
allows us to expect the subsequent experimental obser- 
vation of the novel type of vector solitons and swinging 
dynamics described above. 
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